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AN UNBIASED CODE PHASE ESTIMATOR FOR MITIGATING MULTIPATH 

IN GPS 

FIELD OF INVENTION 



[0001] The invention relates generally to Global Positioning Receivers (GPS) 

receivers or to any navigation receiver that receives pseudorandom noise (PN) encoded 
signals, and particularly, to receivers that that achieve reduction in errors resulting from 
multipath. 



BACKGROUND 



[0002] The Global Positioning System (GPS) was established by the United 

States government, and employs a constellation of 24 or more satellites in well-defined 
orbits at an altitude of approximately 26,500 km. These satellites continually transmit 
microwave L-band radio signals in two frequency bands, centered at 1575.42 MHz and 
1227.6 MHz, denoted as LI and L2 respectively. These signals include timing patterns 
relative to the satellite's onboard precision clock (which is kept synchronized by ground 
stations) as well as a navigation message giving the precise orbital positions of the 
satellites, an ionosphere model and other useful information. GPS receivers process the 
radio signals, computing ranges to the GPS satellites, and by triangulating these ranges, 
the GPS receiver determines its own position and its own internal clock error. 

[0003] The GPS satellites transmit timing information by modulating 

pseudorandom codes on both the LI and L2 frequencies. These codes, referred to as PN, 
or PRN codes, are binary sequences of +1 and -1 as depicted in FIG. 1 (sometimes, l's 
and 0's rather than 1 and -1 are used for notational convenience). Pseudorandom codes, 
known as the P or P(Y) codes, are present on both LI and L2 GPS frequencies; have 
chipping rates of 10.23 MHz; and repeat every week (actually, only the P portion of the 
P(Y) code repeats every week, the Y portion is additional modulation layer for the 
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purpose of anti-spoofing and is unknown to civilian users). The LI channel also has 
modulated upon it a PN code commonly called the C/A code. The C/A code runs at 
1 .023 MHz and repeats every millisecond. The use of PN codes on the GPS broadcast 
signal cause the signal's frequency spectrum to spread or widen, a phenomenon 
commonly called spread spectrum. This makes the signal less susceptible to jamming. It 
also allows GPS satellites to use the same transmission frequencies (LI and L2), but 
different PN codes, a communication technique commonly referred to as code division 
multiple access (CDMA). Most importantly, however, the PN codes provide a means to 
extract timing information from the signal. 

[0004] In acquiring the PN sequence for subsequent tracking, a replica of the 

broadcast PN sequence (See FIG. 2) is produced by the tracking hardware (any PN 
sequence generated by the GPS receiver's tracking hardware, as opposed to the GPS 
satellite, is referred to as a PN model or PN reference). This PN model is then swept 
through time until it aligns with the incoming signal. Alignment is determined by cross- 
correlating the model of the PN sequence with the measured signal. 

[0005] In order to track pseudorandom (PN) sequences in a GPS receiver, a 

delay-lock-loop is often used that employs an early-minus-late (E-L) model of the 
incoming GPS PN sequence. That is, two replicas of the incoming PN sequence are 
produced, whereby one is then shifted l A chip early and the other shifted Vi chip late in 
time. The resulting late code sequence is subtracted from the early sequence. The 
prompt, (the on-time sequence as shown in FIG. 5A), represents the signal that we wish 
to track. A correlation function is defined as depicted in FIG. 3, such that the prompt is 
the reference and is held stationary. When the signal arrives earlier than prompt, its 
resulting correlation function is shifted to the right (delayed) since the signal must be 
delayed before it aligns with prompt. When the signal arrives later than prompt, its 
correlation function is shifted to the left since it must advance to align with prompt. The 
early-minus-late correlation pattern effectively provides the tracking software an error 
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feedback that drives the PN phase generator to track the prompt. This error is known as 
the discriminator of the delay-lock-loop. 

[0006] When the incoming signal arrives early with respect to the model, more 

energy is present in the early correlation result and a positive error is produced. The 
positive error advances the reference phase and the PN (E-L) model moves earlier in 
time. A negative error arises when the incoming signal arrives late and more energy 
shows up in the late sequence. This delays the reference phase. The error continuously 
drives the delay-lock-loop, which thereby adjusts the reference phase generator to 
maintain zero phase error. In a GPS application, the prompt signal is the digitally 
sampled (received) signal corrupted by noise, multipath and filter rounding. Multipath 
signals can cause errors in the prompt and thereby result in a shifted reference. A shifted 
reference causes the delay-lock loop to produce an erroneous code phase. Therefore, 
what is needed in the art is a method to mitigate erroneous effects due to multipath 
signals. 



SUMMARY OF THE INVENTION 

[0007] Disclosed herein in an exemplary embodiment is a method of measuring 

the phase of a pseudorandom (PN) sequence of chips, each chip having a polarity that is 
either positive or negative, and having a time duration, T c hi p . The method comprises: 
generating a reference model exhibiting a reference phase, the reference phase adjustable 
to facilitate alignment with the phase of the pseudorandom sequence; establishing a 
plurality of pulsed-windows over which a plurality of samples of the pseudorandom 
sequence are collected for a selected accumulation interval; and accumulating the 
plurality of samples for each pulsed-window of the plurality of pulsed windows to form a 
plurality of accumulated sums. The method also includes: compensating each 
accumulated sum to form at least one compensated sum, if a number of level transitions 
and non-transitions of the pseudorandom sequence is not equal; and combining the 
compensated sum to determine a phase error from the plurality of compensated sums, the 
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phase error corresponding to a phase difference between the reference phase and the 
phase of the pseudorandom sequence. 

[0008] Disclosed herein in another exemplary embodiment is a method of 

measuring the phase of a pseudorandom (PN) sequence of chips, each chip having a 
polarity that is either positive or negative, and having a time duration, T C hi P . The method 
includes: generating a reference model exhibiting a reference phase, the reference phase 
adjustable to facilitate alignment with the phase of the pseudorandom sequence; 
establishing a pulsed-window over which a plurality of samples of the pseudorandom 

4 

sequence are collected for a selected accumulation interval; and accumulating the 
plurality of samples for each pulsed-window to form at least one accumulated sum. The 
method also includes: combining the at least one accumulated sum to determine a phase 
error, the phase error corresponding to a phase difference between the reference phase 
and the phase of the pseudorandom sequence; and normalizing the phase error by 
dividing the at least one accumulated sum, by a second summation taken only over 
pulsed windows that are based on a number of polarity transitions of the pseudorandom 
sequence. 

[0009] Also disclosed herein in yet another exemplary embodiment is a storage 

medium encoded with a machine-readable computer program code for measuring the 
phase of a pseudorandom (PN) sequence of chips, each chip having a polarity that is 
either positive or negative, and having a time duration, T C hi P , said storage medium 
including instructions for causing controller to implement the abovementioned method. 

[0010] Further disclosed herein in an exemplary embodiment is a computer data 

signal embodied in a carrier wave for measuring the phase of a pseudorandom (PN) 
sequence of chips, each chip having a polarity that is either positive or negative, and 
having a time duration, T C hi P , said data signal comprising code configured to cause a 
controller to implement the abovementioned method. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] FIG. 1 shows a portion of a pseudorandom noise (PN) sequence; 

[0012] FIG. 2 shows the process of aligning a model PN sequence with a 

measured PN sequence. Both sequences are identical other than a time misalignment; 

[0013] FIG. 3 depicts an autocorrelation function that results when two similar 

PN sequences are correlated with one another; 

[0014] FIG. 4 is an example showing the relationship of the phase of a PN 

sequence to the chips of the PN sequence. 

[0015] FIG. 5A presents three PN sequences of similar type and composition, but 

with alignments designated as prompt, early and late; 

[0016] FIG. SB presents the autocorrelation functions of the early and late PN 

sequences with the prompt PN sequence; 

[0017] FIG. 5C depicts the correlation of the difference of a late PN sequence 

from early PN sequences (E-L sequence) with the prompt PN sequence; 

[0018] FIG. 6 shows several chips from a prompt PN sequence, several chips 

from a E-L sequence, and the product of the two sequences; 

[0019] FIG. 7A depicts an E-L properly aligned with the PN transition; 

[0020] FIG 7B depicts an E-L arriving early and the resulting negative error; 
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[0021] FIG. 7C depicts an E-L arriving late and the resulting positive error; 

[0022] FIGS. 8A, 8B, and 8C depict a similar situation as in FIGS. 7A, 7B, and 

7C, respectively, however, with the E-L pulse narrowed to less than one chip width; 

[0023] FIGS. 9A and 9B depict illustrative multipath delayed with respect to the 

direct signal and the combination with the direct signal to form a composite signal; 

[0024] FIGS. 10A and 10B depict how the multipath introduces a bias into the 

track point of the direct signal, causing a non-zero phase error; 

[0025] FIGS. 11A and 11B show how a narrow pulse window (narrow 

correlation) reduces the multipath-induced error; 

[0026] FIGS. 12A and 12B show how a W-discriminator can reduce the multipath 

induced error for multipath of sufficient delay with respect to direct signal; 

[0027] FIG. 13 is an example of the multipath-mitigation windows; 

[0028] FIG. 14 is a block diagram depicting a portion of a GPS receiver involved 

in tracking a PN encoded signal representing in accordance with an exemplary 
embodiment; 

[0029] FIG. 15 is a table indicating total numbers of polarity transitions and non- 

transitions within the 1023 chip C/A code of GPS; 
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[0030] FIGS. 16A and 16B depict illustrative dimensions of a window centered 

near a polarity transition in a PN code for a Low-to-High (L-H) transition and a High-to- 
Low (H-L) transition respectively; 

[0031] FIGS. 17A, 17B, and 17C depict an illustrative PN code, a windowed 

sampling thereof, and sign corrected samples thereof; 

[0032] FIG. 18 depicts sample of points taken over a line segment where spacing 

of samples within the line segment follows a uniform distribution; 

[0033] FIG. 19 shows samples taken on a more realistic PN level transition; 

[0034] FIG. 20 depicts a measured phase error vs. true phase error for one 

embodiment of the present invention; 

[0035] FIG. 21 is a graphical representation of sampling a PN sequence over 

windows centered near the level-transitions and then multiplying the samples by +1 for a 
L-H transition and -1 for a H-L transition; 

[0036] FIG. 22A depicts pulsed windows designed to give phase error curve that 

is less susceptible to multipath induced errors, the pulses aligned with the level 
transition's chip edge as is intended to produce zero phase error; 

[0037] FIG. 22B depicts pulsed windows designed to give phase error curve that 

is less susceptible to multipath induced errors, the pulse advanced to the left sufficiently 
to cause the resulting asymmetric error function to again produce zero phase error; 
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[0038] FIG 22C depicts pulsed windows designed to give phase error curve that is 

less susceptible to multipath induced errors, the pulses delayed to the point that resulting 
phase error is at maximum positive value; 

[0039] FIG. 23 A depicts pulsed windows designed to give phase error curve that 

is less susceptible to multipath induced errors, the pulses advanced to the point of 
minimum resulting phase error; 

[0040] FIG. 23B depicts pulsed windows designed to give phase error curve that 

is less susceptible to multipath induced errors, the pulses delayed to the point of near 
maximum resulting phase error; 

[0041] FIG. 24A depicts the phase error curve that results from using the pulses 

of FIG. 22 and FIG. 23; 

[0042] FIG. 24B provides definition of the dimensions that may be used to 

achieve the desired ratio of areas of the pulses as are need to produce the phase error 
curve in FIG 24A; 

[0043] FIGS. 25A, and 25B depict direct and multipath error signals 

superimposed on one another, for the purpose of showing the multipath reducing benefit 
of the invention; 

[0044] FIG. 26A depicts an improved phase error curve that exhibits better 

immunity to multipath that is delayed by approximately one code chip; 

[0045] FIG. 26B provides definition of the dimensions of pulses used to achieve 

the phase error curve of FIG. 26A; and 
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[0046] FIG. 27 depicts a flow chart illustrating an exemplary process 1 00 for 

phase error computation and utilization of the phase error to control the reference phase 
and to align it with the received PN phase. 

DESCRIPTION OF THE EXEMPLARY EMBODIMENTS 

[0047] A PN code tracker must be able to "reject" the multipath signals to the 

best extent that is possible while still track the desired direct PN signal. The direct signal 
contains the needed timing information that ultimately allows the receiver to compute its 
location 

[0048] Disclosed herein is a new method and system to ascertain a code phase 

error directly, rather than necessarily trying to simply drive it to zero with a track loop. 
Summations are still taken over one or more periodically occurring windows, the size and 
placing of which are tailored to reduce the effects of multipath. Furthermore, a second 
additional summation is taken of absolute values of signal samples over windows 
centered near chip polarity transitions. This second summation is an advantageous 
property of an exemplary embodiment that facilitates normalization. Dividing the first 
set of summations by the additional summation and multiplying by a constant related to 
the window size leads to a measurement of actual phase misalignment between the 
received PN sequence phase and the reference phase. This measurement of phase 
misalignment drives a control loop that keeps the reference phase aligned with the phase 
of the received signal. 

[0049] Second, windows are placed asymmetrically with respect to chip clock 

transitions and such placement results in an asymmetric phase error curve. This 
asymmetric phase error curve is beneficial for reducing the adverse effects of multipath. 
When employing an embodiment of the present invention, the measured-verses-true 
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phase error curve is of unit slope for small displacements of phase error, but becomes 
saturated at a constant positive value to the right (lagging reference phase) and becomes 
zero in a region to the left (leading reference phase). The benefits of the asymmetric 
error function will become more apparent with the disclosure herein. Advantageously, 
the asymmetric error function is less likely to break lock, relative to existing 
discriminators. As a related benefit, the asymmetric error curve is less likely to result in 
inadvertent tracks of a multipath signal rather than the direct PN signal. Tracking 
multipath would mean that the reference phase was aligned closely with a multipath 
signal. 

[0050] Also disclosed herein is a method whereby phase biases are eliminated 

from the asymmetric phase error function by performing additional computations when 
accumulating the phase error to address a disparity in the number of polarity transitions 
(zero-crossings) as non-polarity transitions on the clock edges (integer phase values) of 
the PN sequence. 

[0051] In the discussion that follows, it is assumed that the conventional methods 

of down converting a GPS signal to base-band and conventional digital tracking hardware 
for tracking phase in phase-lock-loop (PLL) are already in place (usually using a phase 
tracking loop called a Costas loop is utilized). Concentration is on PN tracking and not 
the tracking of the signal's carrier. It should be appreciated that carrier tracking occurs 
simultaneously with PN tracking, however, it is well known and not pertinent to the 
understanding of the concepts that of this disclosure. 

[0052] To steer the chip transitions with respect to time the clocking of the model 

PN generator is manipulated. As shown in FIG. 4, a reference PN phase counter 
consisting of two parts: an integer phase 120 and a fractional phase 121, is used for this 
purpose. Both the integer and fractional phase define the total phase of the PN model. A 
delta-phase, or code phase error is usually computed/tracked by a tracking loop and is 



- 10- 



COU-0011 



added to the fractional phase every sample clock period or at some other prescribed rate. 
The continuing accumulation in the fractional phase due to this delta-phase addition 
eventually causes the fractional part to exceed one. At this instant, the integer phase is 
incremented, and an integer's worth of phase removed from the fractional phase. An 
integer phase of zero represents the starting chip of the sequence, and each new integer is 
the next chip in the sequence. The integer count continues until the sequence is 
exhausted (1023 chips in the case of a C/A code) at which point the phase resets back to 
zero. 

[0053] In order to track pseudorandom (PN) sequences in a GPS receiver, a 

delay-lock-loop is often used that employs an early-minus-late (E-L) model of the 
incoming GPS PN sequence. That is, two exact replicas of the incoming PN sequence 
are produced, whereby one is then shifted Vi chip early and the other shifted Vi chip late 
in time. The resulting late code sequence is subtracted from the early sequence. FIG. 5A 
shows the Vi chip early sequence 130, a prompt 131, (or on-time sequence), and a l A chip 
late sequence 132. The prompt, in this case, represents the signal that we wish to track. 
FIG. 5B shows the correlation of the both the early and late sequences with the prompt 
sequence and FIG. 5C shows the correlation of early-minus-late sequence with the 
prompt. The correlation function is defined such that prompt 131 is the reference and is 
held stationary. The axis 134 is such that a delay in time alignment, r , is considered 
positive and to the right. When the signal arrives earlier than prompt, its resulting 
correlation function 136 is shifted to the right since the signal must be delayed before it 
aligns with prompt. When the signal arrives later than prompt, its correlation function 
135 is shifted to the left since it must advance to align with prompt. Note that the early- 
minus-late correlation function 137 crosses zero r = 0 , 133. That is, when early and late 
codes are balanced Vi chip on either side of a prompt code. The early-minus-late 
correlation pattern effectively provides the tracking software an error feedback that drives 
the PN phase generator to track the prompt. This error is known as the discriminator of 
the delay-lock-loop. 
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[0054] When the incoming signal arrives early with respect to the model, more 

energy is present in the early correlation result and a positive error is produced. The 
positive error advances the reference phase and the PN (E-L) model moves earlier in 
time. A negative error arises when the incoming signal arrives late and more energy 
shows up in the late sequence. This delays the reference phase. The error continuously 
drives the delay-lock-loop, which thereby adjusts the reference phase generator to 
maintain zero phase error. In a GPS application, the prompt signal is the digitally 
sampled (received) signal corrupted by noise, multipath and filter rounding. 

[0055] It is instructive to look a little closer at the E-L PN sequence and how it 

* " * * ■* 

correlates with the prompt PN sequence. Consider FIG. 6, which shows several chips 
140 from a prompt PN sequence with each chip having a time-duration of T C hi p . Also 
shown is an E-L sequence, 141, properly aligned in phase with the prompt sequence, and 
the product 142 of the prompt sequence and E-L sequence. Note that although the 
prompt PN sequence takes on values of +1 and— 1, the E-L sequence can take values of 
+2, -2, and 0. However, in practice the E-L sequence is divided by two as it is generated 
so the product of prompt and E-L also takes on values of -H, -1 and 0. 

[0056] Notice that the E-L sequence, when aligned as in FIG. 6, is non-zero at 

chip level transitions of the prompt. Further, the non-zero chips of the E-L are of polarity 
equal to the transition of the prompt. Low-going edges of the prompt produce negative 
E-L chips. High-going edges produce positive chips. Clearly, looking at the product, the 
average value is zero if averaged across time. This is exactly what the hardware 
accumulator does; after forming the product of E-L and prompt, it sums (or integrates) 

the product over a period of time (often 1 to 20 ms for C/A code) called the accumulation 

t 

interval. When E-L and prompt are aligned, the summation is zero. Another way to view 
this is that the E-L chips exactly straddle the chip transitions of the prompt signal Thus, 
the area on the left-hand side of the prompt's transition exactly cancels the area on the 
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right hand side of the prompt's transition. If the E-L signal shifts left (e.g., delayed) or 
shifts right (e.g., advanced), the right-hand and left-hand areas in the product of E-L and 
prompt are unbalanced, and a track error results. In subsequent discussion, the E-L 
sequence will be viewed as a series of "windows". These windows correspond with the 
non-zero portions of the sequence and each window is centered on the prompt's code 
transition when properly aligned and tracking. 

[0057] Using the notation of "windows" centered about the chip transitions, 

consider the correlation of one E-L code window with a signal transition of the incoming 
code. This is shown in FIGS. 7A-7C for three different alignment scenarios: prompt, 
early model, and late model respectively. The top-most diagram in each scenario depicts 
two chips of incoming code (150, 160, 170) with a polarity transition 153 between the 
chips, the diagrams below show the E-L code windows (154, 164, 174) for each scenario, 
and the bottom diagrams show the resulting correlation functions (155, 165, 175). 

[0058] Referring now to FIG. 7 A, the E-L window 154 is centered on the code 

transition 153 and the area on the right 152 balances the area on the left 151. 
Consequently, when the product of incoming code with E-L window is formed, the areas 
sum to zero. The correlation function 155 crosses zero where this occurs 156. The slope 
157 of the correlation function shows how the correlation changes as a result of the 
misalignment. The slope is proportional to the rate at which area is added and subtracted 
due to shifting the E-L code window. The window 164 shifts to the left, as shown in 
FIG. 7B, when the phase of the tracker is advanced in time with respect to the incoming 
signal. In this case, the resulting correlation is negative 166 and this negative error is fed 
back to the tracker. This slows the phase generator of the tracker and causes a delay that 
re-aligns the signals. The opposite occurs as the window 174 is moved right in FIG. 7C. 
Here the correlation is depicted as positive 176. 
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[0059] Ideally, the code tracker is configured to keep the window centered and 

consequently, only the slope of the correlation function near the zero crossing is of 
significance while maintaining track. This means that once a stable track is obtained, the 
window may be narrowed, if desired, as depicted in FIGS. 8A-8C. The slope 202 of the 
correlation functions on the bottom row of FIGS. 8A-8C is the same as those of FIGS. 
7A-7C, but in FIG. 8A-8C we see that the correlation function 201 flattens on the top and 
bottom. The slope is the same because the narrower window "picks up area" at the same 
rate as the wider window as it is slid across the edge of the code-transition. 

[0060] Narrowing the window as described has commonly become known as 

"Narrow Correlation". One advantage of narrow correlation is that it improves the Signal 
to Noise Ratio (SNR) of tracking. It will be appreciated that in an actual receiver, the 
incoming signal may be corrupted with noise. In tracking, it is the edge transitions that 
are important. By narrowing the amount of data gathered away from the edges, we 
essentially eliminate data that is not providing additional information, but which is 
adding noise to the tracking discriminator. 

[0061] The problem with E-L discriminators, and even the Narrow Correlation 

E-L discriminators, is that they are affected by reflected signals from nearby objects, 
these reflections known as multipath (i.e., the signal travels multiple paths to get to the 
receiver, the direct path and via the paths that result from reflections off objects nearby). 
When multipath is present, there is more than one PN signal arriving: the desired direct 
signal and one or more reflected signals. The reflected signals always arrive late with 
respect to the direct signal. 

[0062] The reflected signals are typically weaker than the direct signal. 

Furthermore, the arriving reflected signals are sometimes reversed in polarity to the direct 
signal. This is due to: 1) the reflection flips the phase 180 degrees, and 2) the phase of 
the carrier varies with path distance. For every 19 cm of path difference, the wavelength 
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of the Li's carrier, the phase will rotate through 360 degrees when tracking C/A code for 
example. The multipath can arrive at any phase between 0° and 360° with respect to the 
direct signal. The "I" and "Q" channels that feed a typical Costas carrier-tracking-loop 
will see the polarity of the multipath PN sequence flipped whenever its carrier phase is 
more than 90 degrees offset from the direct signal 

[0063] FIG. 9A shows an example of both a direct 210 and a reflected (multipath) 

signal. In this figure, the reflected signal 211 exhibits the same polarity as the direct 
signal but is delayed by 3/8 of a chip and exhibits 3/8 of the amplitude as well. The 
combined direct and reflected signal 212 is shown at the bottom of the figure. FIG. 9B 
shows a similar situation, however, the reflected signal 214 is opposite polarity to the 
direct signal. For clarity, we show only one reflected signal in each figure, however, in 
practice, there are often multiple reflected signals arriving at varying delays, amplitudes, 
and polarities. The discussion presented herein easily extends to multiple reflected 
signals. 

[0064] Now consider the impact of a signal that contains multipath on an E-L 

discriminator. As shown in FIG. 10, the E-L window 223 must still balance the area on 
the right and left hand sides. However, this causes the E-L window 223 to be delayed 
with respect to the direct signal when the multipath 221 is the same polarity as the direct 
signal. The E-L window, 227, advances from the direct signal when multipath 225 is of 
opposite polarity to the direct signal. 

[0065] What can readily be observed is that due to the shape of the combined 

signal, the effect of multipath that is opposite polarity to the direct signal is more 
pronounced than the effect of multipath of the same polarity as the direct signal. That is, 
the advance caused by reversed-polarity multipath is larger than the delay caused by 
multipath of the same polarity as the direct signal, assuming the multipath is identical 
except for its polarity. r 
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[0066] Narrow correlation provides some performance improvement. Consider 

Figures 1 1 A and 1 IB, in this instance, the largest error caused by the multipath cannot 
exceed Vi the window width of the narrow correlator. Narrow correlation is impacted by 
multipath and is not much better than a conventional E-L correlator for multipath that is 
from nearby reflection points. 

[0067] The utilization of a narrow window may readily be taken one step further. 

In one instance, a different shaped window, known as a W-discriminator is discussed. 
Employing a discriminator with this window shape will eliminate the error introduced by 
multipath completely, provided the time-delay of the multipath signal is sufficiently large 
that the code transitions in the multipath lie beyond the boundary of the window. FIG. 12 
shows a W-discriminator window waveform 246. Note that the window is essentially an 
even symmetric narrow correlation type window with two side lobes of opposite polarity 
to the main window at Vi the amplitude. The side lobes 251 and 252 are situated on either 
side of the main window 250. 

[0068] It may readily be appreciated that rather than trying to analyze the 

combined signal 245, it is beneficial to analyze the direct signal 240 and reflected 
multipath signal 244 separately. Clearly, the W-discriminator' s area 243 averages to zero 
and, as such, its product with a constant waveform will average to zero as well. When 
the multipath 244 is sufficiently delayed, as it is depicted in FIG. 12, the multipath 
appears constant over the extent of the W-discriminator window. Consequently, the 
multipath has no effect on the average value of the correlation. Analyzing the direct 
signal 240 it may be seen that the lobes of the window are such that the product of the 
lobes with the direct signal has area 241 and area 242 of equal size and opposite sign. 
This holds true for shifts of up to Vi the width of the main window 250 of the W- 
discriminator. So there is no impact of the side lobes on the direct signal when the W- 
discriminator window is nearly centered with the code transitions. This leaves the main 
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window, which has an effect similar to narrow correlation on the direct signal. That is, 
shifts to the left or right cause increasing or decreasing area in the correlation. Therefore, 
the W-discriminator is a narrow correlation type discriminator that is immune to a 
sufficiently delayed multipath. 

[0069] Therefore, it may also be appreciated that the drawbacks of the W- 

discriminator include: First, it is easier to lose lock, since the W-discriminator pattern 
goes to zero after relatively narrow windows on either side of the discriminator function. 
Furthermore, this zeroing of the discriminator makes it possible that the control loop 
would zero out the direct sequence and thereby, track the multipath. Second, the W- 

r 

discriminator includes added side lobes that are not present, in a narrow correlator, the 
additional window area results in more noise being picked up in the correlation process. 
Since all information needed for accurate tracking is in the code transitions of the 
incoming PN signal, any extra regions that are added to the correlator introduce noise, but 
no additional information. 

[0070] An even more generalized version of the windows needed for tracking 

with reduced multipath are known as Multipath Mitigation Windows (MMWs). They 
generalize to include the W-discriminator window, but take into account additional 
window shapes. Here, asymmetric windows are introduced that have the potential to be 
narrower than the blanked PRN window since they are less likely to lose lock. 

[0071] The MMW windows are formed from combinations of the original PRN 

code, where these combinations are shifted in time, scaled in sign and amplitude, and 
then only regions near the edge transitions are utilized. Due to the fact that the windows 
are sometimes asymmetric in nature, they no longer balance across a single chip 
transition and you must consider the entire sequence of chips over which the correlation 
takes place. This windowing strategy assumes that there are roughly as many transitions 
as non-transitions in the C/A code signal and that asymmetric windows can balance 
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transitions with non-transitions. This results in a MMW that are shaped like FIG. 13, for 
example. 

[0072] Unfortunately, asymmetric windows may introduce a tracking bias 

because there is not an exact balance of transitions with non-transitions. Transitions in 
the C/A code can differ from non-transitions by as much as 65 (see FIG. 15) out of 
sequence length of 1023 chips. Moreover, asymmetric MMW windows are placed at all 
clock edges of the code, not just where the transitions take place. This adds noise and 
may degrade tracking performance. 

|0073] Referring to FIG. 14, a block diagram is shown representing a portion of a 

GPS receiver 10 in accordance with an exemplary embodiment. The PN encoded signals 
from a plurality of GPS satellites are received by an antenna 300, and then supplied to an 
RF Front End and Down Converter. Down converted data is digitally sampled, the 
sampled data 302 of which is supplied to multiple GPS processing channels. FIG. 14 
depicts one such channel, all of which perform a similar function. For simplicity, 
additional paths of processing that are sometimes included in advanced GPS receivers are 
omitted from FIG. 14. The block diagram, as shown, pertains to processing of the LI 
signal component of GPS, although the inventions disclosed herein readily apply to both 
LI an L2 GPS signals and even other Navigation satellites such as the proposed Galileo 
system. For each channel, an I-reference generator 303 generates an "in-phase" signal 
that is the same frequency and phase as the carrier component of the sampled data. A Q- 
reference generator 304 produces a "quadrature-phase" signal that is ninety degrees out- 
of-phase to the carrier component. Both the I and the Q reference signals multiply the 
incoming data samples for the purpose of removing the carrier component. In an 
exemplary embodiment, focus in on the I-channel 305 by way of illustration since when 
the receiver is phase-locked, most of the signal energy is in the I-channel. It should be 
appreciated of course, that this is not considered a limitation, however. The embodiments 
disclosed herein are equally applicable to the Q-channel, and to combinations of I- and 
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Q-channels. The Q-channel, 306, is shown for completeness. The Q-channel 
accumulates and dumps periodically a Q-prompt signal 307. This Q-prompt signal is 
processed by a Phase-Lock-Loop (PLL) 310 that drives a Carrier DCO 311 to maintain 
phase-lock on the carrier of the incoming data stream. Also shown is a Q-track signal 
308 (shown un-connected) that is used primarily for control loop "pull-in" prior to 
obtaining phase lock. In addition, FIG. 14 depicts an I-prompt signal 309 that aids the 
PLL by providing the proper sign to the Q-prompt signal, a step that is employed to 
address half-cycle phase ambiguity. One such PLL that operates with a half-cycle phase 
ambiguity is commonly known as a Costas loop. The receiver 10, in subsequent 
processing of the I-prompt signal, extracts navigation data that has been modulated on the 
GPS broadcast signal and that contains relevant information such as the location of the 
GPS satellites. 

[0074] The I-channel includes, but is not limited to, a PN encoded signal that is to 

be tracked. In an exemplary embodiment, tracking this PN encoded signal is 
accomplished by measuring a phase alignment error between the phase of the received 
PN signal and a locally generated reference phase 324 utilizing an approach that yields a 
true measure of the tracking error. An E-L code generator, 312, and a prompt code 
generator, 313, are employed, but serving different purposes from their conventional 
application. In an exemplary embodiment, the Prompt and E-L codes from these code 
generators are supplied to pulsed-signal generators within block 314. This aids in the 
purpose of selecting proper pulse windows from a set of pulse generators 315, 316 and 
further facilitates the determination of the polarity of such pulses. It will be appreciated 
that although two pulse generators are shown for illustration of an exemplary 
embodiment, other implementations including fewer or a greater number of pulse 
generators is possible. The pulsed signals are then multiplied by the I-channel data and 
fed to a set of accumulators, 318, and 319. The accumulators provide running sums of 
the pulse-multiplied data, which is dumped periodically (for example every millisecond) 
to another processing block 321. The processing block 321 combines the accumulated 
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data and normalizes it to produce a code phase error 322. One feature of the 
normalization process disclosed herein 321 is that the accumulated data is scaled to 
remove biases that would otherwise occur in the code phase error 322 due to the 
asymmetric approach taken when forming this error. The code phase error 322 is 
supplied to a Code Track Loop, 323, the purpose of which is to provide delta-code values 
325 that are summed within a code DCO 326. The end result is that the reference code 
phase 324 is kept synchronous with the code phase of a PN sequence received from a 
GPS satellite and which is present in the I-channel 305. The reference code phase 324 is 
required when forming a pseudo-range, as is commonly known in the art, and this 
pseudo-range, when processed together with pseudo-ranges from the plurality of tracked 
GPS satellites, allows the location of the receiver to be determined. Furthermore, 
achieving synchronous reference code phase facilitates alignment of a PN model with the 
actual PN signal. Such alignment is also utilized for further, well-known aspects of GPS 

f 

data processing including, but not limited to, data-demodulation, carrier phase tracking, 
and signal strength monitoring. 

[0075] The reference code phase 324 is steered by a conventional feedback 

control loop 323 in such a manner that it remains aligned with the received PN code 
phase. The control loop is driven by a code phase error term 322. The code phase error 
term 322 is an estimate of the phase difference between the reference code phase and the 
received code phase and is defined to be positive when the reference code phase lags the 
actual code phase, zero when the phases are aligned, and negative otherwise. This 
convention for the sign of the code phase error 322 is consistent with standard feedback 
control practice, but the sign is not important, as long as the control law is defined such 
that, ultimately, an observed lag in model phase results in increasing the rate of phase 
accumulation so that the reference code phase 324 "catches up" to the incoming actual 
code phase. Conversely, the reverse should occur when the reference code phase 324 is 
leading the incoming signal phase and the control law should slow the accumulation of 
reference code phase 324. 
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[0076] In an exemplary embodiment an unbiased, direct measure of code phase 

error 322 is produced that is not corrupted by multipath generated from objects in the 
surrounding area of the antenna 300, provided the signals are delayed sufficiently with 
respect to the direct signal. Such multipath can cause detrimental effects in the position 
accuracy of a GPS receiver 10. 

[0077] It should be noted that to determine the phase of the incoming PN signal, 

it is necessary to measure the occurrence of some feature within the signal that provides a 
distinguishable mark of the signal's progression through time. It is the transitions in 
polarity of the PN signal that provide such a mark. Such transitions occur only on a 
whole integer value of the phase of the PN sequence. As stated earlier, the phase of the 
sequence of PN chips is a number having both a whole integer portion and a fractional 
portion. The whole integer portion is a count of completed chips of the PN sequence, and 
the fractional portion is the fraction of the chip that is in progress, either partly 
transmitted or received. At the beginning of each new chip, the phase, by definition, is 
exactly an integer value (i.e., the fractional part is zero). Once again, FIG. 4 provides an 
example depiction of the phase of a PN sequence after 2.5 chips of the sequence has 

j 

transpired. In the example, the integer part 120 is two chips in duration and the fractional 
part 121, is 0.5 chips. 

[0078] Locally, the receiver 10 generates a reference code phase 324 to follow the 

true code phase of the received sequence. If the reference code phase 324 is aligned with 
true code phase, it should achieve an integer value simultaneously with a chip transition 
of received sequence. If not aligned, there is an error in the reference code phase 324 as 
compared to the code phase of the received signal. As indicated, the error is used in a 
control loop to correct the reference code phase 324. The difference in time between 
incoming polarity transitions and reference phase integer values, when divided by the 
chip length, is the desired code phase error 322. 
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[0079] It is assumed that the reference code phase 324 is initially aligned to the 

incoming PN sequence by commonly known GPS search practices and that it is 
maintained by subsequent control-law feedback such that the phase error is a fraction less 
than one chip, and is usually significantly less than one chip. Therefore, we desire only a 
fractional phase error. 

[0080] As a PN code is clocked out of a PN code generator, there are exactly four 

events that may transpire with respect to the polarity of the code across a clock boundary. 
The code can either change state or not, and the new polarity can be either +1 or -1 . The 
four possibilities are defined as: 

H-H Non-Transition - The code is +1 prior to the clocking of a new chip, and +1 
afterwards; 

L-L Non-Transition - The code is -1 prior to the clocking of a new chip, and -1 
afterwards; 

H-L Transition - The code is +1 prior to the clocking of a new chip, and -1 
afterwards; and 

L-H Transition - The code is -1 prior to the clocking of a new chip, and +1 
afterwards. 

[0081] The table depicted in FIG. 15 provides the totals of the various clock-edge 

events that occur within the C/A code of GPS. 

[0082] It should now be appreciated that there are only two possibilities for level 

transitions, the H-L and the L-H transition. FIGs. 16A and 16B depict a detailed view of 
one L-H transition one H-L transition for illustration. In FIG. 16A, the point in time 403 
at which the reference phase becomes an integer and at which a new chip would be 
clocked is labeled *C\ The nomenclature 'C is used to indicate a centerline of the 
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reference model's PN level transition. A window centered on 'C and extending a 
distance of 4 w' on each side of the centerline 'C is shown, having left edge 401 and right 
edge 402. The edge of the polarity transition 408 of the sequence 405 occurs prior to the 
centerline 'C so that the reference model lags the incoming signal by code phase delay 
'd\ Although the transmitted sequence originates as a series of chips with levels of +1 
and -1, the received signal is shown to have levels of +A and -A (407) where the 
amplitude ' A' is typically different from one. This is due to a variety of factors, 
including broadcast power level, signal attenuation, hardware processing gain etc. There 
is actually no loss of generality when assuming a signal level of +1 and -1, but we use 
+A and -A in FIGs. 16A and 16B to depict a more general and the real situation. 

[0083] A code phase discriminator is developed based upon the area enclosed by 

the window borders 401 and 402, and the signal 405 of Fig 16A. The convention shall be 
taken that a negative signal level produces a negative area and a positive signal produces 
a positive area. The total area bounded on top and bottom by the signal 405 and 
contained within the window shown spanning from 401 to 402 is thus: 

total area = = A(w + d) - A(w - d) = 2 Ad (1). 

[0084] Now consider one more computation of the area, however, this time using 

the absolute value of the signal. If the absolute valued of the signal is taken when 
computing areas, we shall refer to the area as the "absolute area". 

[0085] In FIG, 1 6A for instance: 

absolute area = = A(w + d) + A(w - d) = 2 Aw (2). 

[0086] Now, dividing the total area by the absolute area: 
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total area / absolute area = 2Ad I 2Aw = d/w (3). 

[0087] The desired delay, d, is then: 

d = w (total area / absolute area) (for L-H transitions) (4). 

[0088] Performing the same analysis on FIG. 16B, we arrive at a similar result, 

except for the presence of a minus sign: 

d = -w (total area / absolute area) (for H-L transitions) (5). 

[0089] It may readily be appreciated then that the code phase delay, d, may be 

computed by performing area sums and ratios about either the H-L transition or the L-H 
transition. In FIG. 16, a positive code phase delay, d, was depicted where d has been 
defined positive whenever the reference code phase lags the phase of the received 
sequence. When the reference code phase e.g., 324 leads the phase of the received 
sequence, d is negative, 

[0090] Although, so far, the waveforms presented have been shown as continuous 

functions of time, the actual processes for these waveforms is digital (implemented in 
DSP, ASIC, FPGA, discrete logic or other suitable digital hardware). The hardware 
samples the analog signal at a constant rate but is asynchronous to the code and carrier 
transitions. FIG. 1 7 depicts the continuous and sampled waveforms for several chips of a 
PN code. FIG. 17A shows a continuous representation of the incoming PN code, denoted 
S(t). FIG. 17B shows the samples 417 of S(t) with these samples denoted S(k). As is 
common practice, the argument in parenthesis, denotes that S(t) is a continuous 
function of time. The argument indicates that S(k) is sampled data taken at discrete 
points in time and the value of is an index associated with the discrete point. Vertical 
tick marks 414, on the horizontal axes depict the sample times. Samples occur at fixed 



-24- 



cov-oou 



intervals denoted T s , 415. Samples S(k) are taken only in windows 410 that span a region 
that encompasses the chip edge. The window center may or may not coincide with the 
actual chip edge. In FIG. 17 for example, the windows are depicted slightly off to the 
right of center of the actual chip edge. Due to the window being slightly to the right, 
more samples arise following (to the right of) the actual chip edge than prior to the chip 
edge. 

[0091] R i s desirable to normalize the delay, d, and remove its dependence on the 

transition type, H-L or L-H, as indicated in the change of sign between Eq. (4) and Eq. 
(5). FIG. 17C shows another signal, denoted S\k) , which is produced by performing the 

following logic. For each H-L transition 412, reverse the polarity of S(k) by multiplying 
by -1 and for each L-H transition 413 keep the existing polarity. This results in each 
window of data being multiplied by either +1 or -1 as shown by 418 and 419. 

[0092] The sample period 415 is chosen preferably to not evenly divide the chip 

period 420 of the PN code. Furthermore, it is chosen not to divide evenly any N chips 
where N is any number ranging from one to some relative large number (usually more 
than several hundred). This is assured for the PN code's nominal frequency, and also for 
slight perturbations from the nominal frequency due to Doppler shift arising from either 
satellite motion, receiver motion or receiver clock bias. When the sample period is 
chosen in this manner, the samples S(k) fall at different points within each of the 
windows 410. That is, the first point in each window, and therefore each subsequent 
point, will be offset relative to the edge of the window by an offset that is unique to each 
different window. The occurrence of samples, relative to the edge of the window, is thus 
treated as a statistical outcome of a uniformly distributed process. 

[0093] In one implementation of an exemplary embodiment, at least 512 windows 

worth of data are processed to generate a single estimate of code phase error. The large 
number of windows assures that properties relating to the statistical nature of the uniform 
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distribution assumption can be applied in a useful manner. This is exemplified by FIG. 
18. Here, a line that extends from 430 to 432 represents a window of time of duration W 
{Wis the length of the segment). Points such as 433, 434 and 435 along the line represent 
sample points taken according to a uniform distribution of such points along the line. By 
the laws of statistics, on average, a certain fraction of the total points that lie within the 
entire window of time 430 to 432 will also fall within the window segment of length x 
spanning 430 to 431, and this fraction will be equal to the ratio x/W. The terminology 
"on average" makes clear that this property does not necessary hold true for any one 
window's worth of samples. But as more and more windows are taken together and 
averaged, the ratio of the number of samples falling within the region 430 to 431 to the 
number of samples falling within the total window 430 to 432 becomes increasingly 
likely to fall close to x/W. Sampling over multiple windows produces a statistical 
outcome that is similar to averaging. Thus, considering an entire ensemble of windows, a 
count taken of all samples that fall within front portions 430 to 431 of the windows can 
be divided by a count of all samples within all windows to achieve a ratio approaching 
x/W. It will be appreciated that the certainty in the ratio being near x/W increases with 
increasing numbers of windows (or samples). 

[0094] Thus far, the emphasis of the disclosure presented has been on the ratio of 

number of samples in one segment of a window to the number of samples in the entire 
window and it has been noted how the confidence in this ratio improves by summing 
over multiple windows. But the abovementioned property readily extends to the values 
of the samples themselves if areas, rather than line segments are used to define the ratio. 
That is, under certain assumptions, the ratio between two summations of samples taken 
from different portions of a series of windows approaches the ratio of areas traced out by 
the continuous signal S(t) in each portion of the window. Moreover, the ratio is 
increasingly refined as the number of windows and thus the number of samples increases. 
The assumption is that the data, over which each of the two summations is taken, follows 
similar trends in each window as shown in FIG. 19. Here, the windows 440 and 441 are 
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centered on similar regions, specifically the L-H transitions. By virtue of the fact that 
sample values 444 and 445 and others follow a characteristic shape of the waveform 442 
and 443 of the PN code as it makes the L-H transition, they are roughly consistent from 
one window to the next. Furthermore, the uniform distribution of sample points along 
with a large number of samples assures that samples are evenly and densely dispersed 
across the signal S(t). Note that the data waveforms 442 and 443 have been shown with 
rounding indicating the effect of filtering in the RF front end of the receiver 10. Even 
though the waveforms 442 and 443 are not identical, they are similar enough in shape and 
amplitude for the arguments of consistency to hold true. The values such as 444 and 445 
of the waveform are consistently of similar size and follow similar trends within each 
window and the densely dispersed samples, on average, characterize this consistency. 
Thus, as stated, the ratios of the windowed sums of sample values tend to equal the ratios 
of their corresponding areas under the S(t) curve. 

[0095] As a further point, it will be appreciated that there are benefits of having 

multiple windows in the accumulation interval. However, it may also be appreciated that 
only one window may be utilized rather than a series of windows, provided that enough 
data points were taken within the single window to accurately fit the curve S(t). In 
practice, however, the data is corrupted by noise. Furthermore, this noise is highly 
correlated for points that are near to one another. By sampling over multiple windows 
that are spaced relatively far apart (e.g., one chip length) the noise is substantially de- 
correlated. Once noise is no longer correlated, the inclusion of additional points in a 
summation has the effect of averaging out and thereby canceling the noise. That is, the 
effect of the noise is removed to a high degree so that the summation more accurately 
reflects the true properties of the noise-free PN waveform over the window. 

[0096] Returning now to Eq. (4) and Eq. (5) that were derived using areas as 

depicted in FIG. 16. in light of the observation that summations are proportional to areas, 
summations may be substituted for areas. That is: 
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d - wx 




for samples taken over L-H transitions 



(6); 



d = - 




for samples taken over H-L transitions 



(7); 



where samples e LH windows denotes that samples are taken over the L-H windows and 
samples e //£ Wovvs denotes that samples are taken over the H-L windows. The actual 

number of H-L or L-H windows over which summations occur depends on a pre-selected 
accumulation interval. In a C/A code tracker of an exemplary embodiment, the number 
of windows that fall within the C/A code sequence of 1023 chips that defines the 
accumulation interval. 

[0097] As a further processing step, all data from the H-L transition is multiplied 

by (-1), and then summed over both types of transitions, treating the summation as one. 
Advantageously, the S\k) series as depicted in FIG. 17C facilitates this computation. 

Letting S(k) denote the original samples within the H-L and L-H windows, then S f (k) is 

derived from S(k) by multiplying any data from S(k) that falls within a H-L transition by 



(-1). 



That is: 



S'(k) = S(k) if k e LH 



window ' 



: and 



S\k) = -S(k) tfkeHL 



window ' 
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The code phase error, d, is then given by: 

d = wx keLH !^^ , (8). 

where the samples A: are taken from both the L-H windows and H-L windows as indicated 
by k e LHHL yvindows . It should be noted that these windows are of width D = 2xw. The 

windows are separated from one another by an integer multiple of one code chip's worth 
of phase when the exact multiple depending on the occurrence of L-H and H-L transitions 
within the particular PN code's waveform. 

[0098J Therefore, disclosed herein in an exemplary embodiment is a way to 

directly measure the code phase error d, 322 (Figure 14) by utilizing summations of data 
taken from windows roughly centered on transitions of the PN signal. The code phase 
error d, 322 when supplied to a feedback control law, provides fine-tuning of the 
centering of the window by causing adjustments to the reference model phase 
accumulation rate. Note that the units of d are also the units of w. If w is expressed in 
chips (e.g., one whole phase increment) then the code phase error d is expressed in chips 
of code phase error. If w is in nanoseconds, then so is d. 

[0099] One of the advantages and benefits of the disclosed method is that, in an 

exemplary embodiment, the computed code phase error, d as computed herein is not 
simply proportional to the code phase error (as in existing delay-lock-loops of current 
systems), rather d is an actual determined value for the code phase error, 322. 
Furthermore, d is a direct measure of code phase error regardless of signal strength of the 
received PN code. This means that the control law for maintaining code lock is easier to 
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design and implement since there is not a dependency of loop gain on the received signal 
strength. 

[00100] It will be appreciated that Equation (8) is accurate as long as the code 
phase error, d 322 does not exceed the half window width, w. Once the code phase 
exceeds w, Eq. (8) will no longer produce an accurate measure of code phase error d 322, 
but will instead saturate at +w or -w. This characteristic is depicted in FIG. 20 where 
measured code phase error, d 322, is plotted versus true phase error. The plot 450 is of 
unit slope (i.e. a one to one ratio) for true phase errors less than w (451) in magnitude and 
remains constant for true phase errors in excess of w. Nevertheless, the sign of the 
measured code phase error agrees with the sign of the true error and this holds true for up 
to one chip's worth of error. Because the sign is correct, a feedback control law 
employing this error works correctly to null out the error. With proper sign, the feedback 
control adjusts the model reference phase to the point that the window eventually 
becomes centered on the transition where the phase error measurement is reliable. 

[00101] It will be appreciated that to accommodate larger error, the window could 
be widened, but the phase error measurement would then be increasingly affected by 
multipath as the window widens. 

[00102] Consider multipath received signals that are sufficiently delayed with 
respect to the direct signal that its polarity transitions are outside the window. The affect 
of this multipath is to alter the result of the summations in Eq. (8). Due to the delay of 
the multipath, on average, the multipath signal maintains constant amplitude throughout 
the window. Let m(k) denote the sampled multipath signal and let M be its average 
amplitude. Furthermore, let N denote the total number of points within the summations 
of Eq. (8). Assume that the window is exactly centered on the level-transition as is often 
the desired situation. This means that the term ^ S\k) should be zero. If the average 

amplitude of S(k) is A, then Eq. (8) reduces to: 
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multipath 



wx y m(k) 

keLHHL^ nd0W 



^\S(k) + m(k)\ 



wx N x M wxM 



Nx(A + M) A + M 



if the multipath has the 



same polarity as S(k); and 



multipath 



wx y m(k) ■ 

keLHHL^ ind0WS 



Z\S(k)-m(k)\ 

kf~LHHL^ indow 



wxNxM wxM 



Nx(A-M) A-M 



if the multipath has the 



opposite polarity to S(k); 



where d mu iti P ath is the erroneous code phase error measurement that results from the 
multipath. If M is fairly small compared to A, then the above equations both simplify to: 



^ multipath 



wxM „ w 
forM«A 



[00103] It may be seen that multipath induced error depends on the ratio of the 

multipath amplitude M to the signal amplitude A. Furthermore, the error depends directly 
on the window's size w. Advantageously, decreasing window size decreases the effect of 
multipath on the code phase error measurement, however, the window size must still be 
maintained large enough to avoid having the code tracking loop lose lock. A very small 
window substantially reduces the track loop gain as track errors become large. 



[00104] Returning to Eq. (8), FIG. 21 illustrates a visualization of the forming of 

S\k) graphically. Recall that S'(k) is defined as the data samples that fall within 

windows of the received PN code, with each window centered, or nearly centered, on a 
level-transition. Furthermore, data that is sampled in the vicinity of the H-L transitions 
are multiplied by -1 to reverse their polarity when forming the samples S'(k) . The 

received PN sequence is shown as 460. Samples of this sequence are taken at fixed 
intervals represented by the tic marks such as 461. Data is sampled only in windows 
such as 462 and 463. Furthermore, data sampled in windows surrounding a H-L 
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transition is multiplied by -1. The window 462 surrounds an H-L transition and is shown 
extending downward to -1 . The window 463 surrounds a L-H transition and extends 
upward to +1 . To make the distinction that these windows have an associated polarity, 
we refer to them as "pulses". If the composite of all such pulses is used to form a new 
signal 464, it will take values of +1, -1 and zero. The signal 464 can be multiplied by the 
samples of the signal 460 to produce the series £" (k) . The signal 465 outlines the shape 

of S\k) . Note that when the pulses are exactly centered on the chip transition 467, the 
code phase error denoted in this figure as 466 reduces to zero. At this point, the 
summation of samples of S\k) , if taken over a suitably long accumulation interval, will 
be zero since the negative and positive portions will cancel. 

[00105] Now consider a means to reduce or eliminate multipath in accordance with 

an exemplary embodiment. Let another signal u\j)be introduce and summed into the 

numerator of Eq. (8) such that »'(./)] = 0. This leads to: 

^ k^LHHLnindnvs J e A 234 wdow /Q\ 

keLHHL uindows 

[00106] Since [JV (j)] = 0, Eq. (9) yields the same result as Eq. (8). The data 

series u\j)is generated in means similar to the generation of S\k) . To form u\j) , the 
received PN code is multiplied by a set of windows that we denote the A234 windows . Like 
the LHHL^ nd0WS , the A234 windows have an associated set of pulses. When these pulses are 
multiplied by the received PN sequence, they produce a balanced result that sums to zero. 

[00107] FIG. 22 depicts additional pulses (i.e., areas associated therewith) added 

for the purpose of reducing or eliminating multipath. The composite pulsed signal 501 
depicts a set of pulsed- windows, which have been shown extending positive or negative 
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to represent their polarity. The dotted line, 500, represents the received PN code 
sequence. Both a non-transition at clock edge 506, and a L-H transition at clock edge 
507 are shown for the received PN code. A H-L transition is not shown, with the 
understanding that signs of the pulsed-windows reverse in the case of a H-L transition. 
The pulse 505 of absolute area, Ai is the previously discussed LHHL window. Additional 
pulses 502, 503, and 504 of absolute areas A3, A2 and A 4 , respectively have been added 
to form the composite signal 501. The pulsed window 502 of area A3 directly precedes a 
non-transition portion of the PN code as depicted, and always takes the same sign as the 
reference PN code. Consequently, the product of A3 with the reference PN code is 
always positive. If the reference code phase is aligned with the received PN code's phase 
to within roughly one chip, the product of A3 with the received PN code is also positive 
as indicated by the sign in 502. It remains positive for up to a chip width of 
misalignment between the reference model signal and the received PN code. 

[00108] Products of pulses Ai, A2 and A4 with the received PN code can take 

either positive or negative values according to the alignment of the composite pulsed 
signal 501 with the received PN code 500. In FIG. 22A, the composite pulsed signal is 
properly aligned with the PN code to achieve zero code phase error. By design, it is 
aligned as in FIG. 22 A when "perfect" phase tracking occurs and is designated as being 
"on time" since it is neither early (advanced) nor late (delayed). In this situation, the 
product of the pulses with the received PN code results in signs shown in 502, 503 and 
504, which are either *+' or During "on time" alignment, the product of the pulse Ai 
with the PN code averages to zero since it is balanced on either side of the code 
transition. Consequently, no sign is given for 505 in FIG. 22A. 

[00109] Turning to FIG. 22B, the composite pulsed signal 508 for the reference 
arrives early with respect to the received PN code and has been sufficiently advanced that 
all pulses lie to the left of the L-H transition of the received code. The signs that result 
when the product is formed are then positive for A3 and A4 and negative for Ai and A2. 
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In FIG. 22C the composite pulse signal 509 arrives late with respect to the received PN 
code and has been sufficiently delayed so that all pulses lie to the right of the L-H 
transition of the received code. Products of Ai and A2 with the received PN code are now 
positive and A4 becomes negative. As is always the case for small phase offsets, the A3 
product remains positive. 

[00110] In an exemplary embodiment, Ai and A2 are of selectable areas, each area 

defined by its respective width. The introduction of area A2 causes an asymmetric effect 
that may be seen in the phase error graph of FIG. 24A. It will be appreciated that this 
asymmetric nature of the phase error graph is advantageous for reducing multipath while 
maximizing ability to maintain lock. 

[00111] Areas A3 and A4 are related to A\ and A2 by the following 

A 3 =A 2 + A } /2 
A 4 =A,/2 

Let 2 represent the sum area of the product of the composite set of pulses with the 
received PN sequences. For the three cases depicted in FIGs. 22A - 22C: 

(a) On time: Z = A 3 - A 2 - A 4 = (A 2 + A l 1 2) - A 2 - A { 1 2 = 0 ; 

(b) Advanced: 2 = 4 + A 4 - ^ - \ = + A l 12) + 12 - 4 - 4 = 0 ; and 

(c) Delayed: Z = A 3 + A x + A l - A 4 = + A x 12) + A x + - A l 12 = A l + 2A 2 . 

[00112] It should be noted that for both the "on time" and the "advanced" cases, 

the sum area is zero. The "delayed" case yields a positive result. It should be 
appreciated that the placement and shapes of the pulses of areas Ai, A2, A3 and A4 of 
FIG. 22 are not arbitrary. Their placement has been chosen to yield desired 
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characteristics in the resulting phase error curve as depicted in FIG. 24A. Furthermore, 
consideration has been given to the fact that in a given receiver 10 the processing of the 
PN signal disclosed herein is "down stream" of filtering on the PN signal that causes 
rounding of the PN code edges as well as transitions of polarity that are not 
instantaneous, but sloping. When the pulses are properly phase aligned with the 
incoming PN sequence, the pulse Ai cancels areas on either side of a sloping transition as 
well as an abrupt transition. Furthermore, pulses A2 and A4 are far enough removed from 
the transition, that they should see similar signal levels as seen by pulse A3 and the three 
pulses should cancel when multiplied by the data. 

[00113] In view of the need for computing the true (e.g., deterministic) code phase 
error rather than a proportional phase error, it is beneficial to divide all sum areas by the 
absolute value of the areas Ai + A2, as well as multiply by the product ofD/2 where D is 
the combined width of Ai + A2. That is: 



accurately be measured are depicted in FIGS. 23 A and 23B. Here the maximum shifts of 
the pulse signal are shown for which the measured phase error remains accurate and 
linear. From the figure, it may be observed that: 

(a) Advanced to Edge of Linear Region (minimum measurable phase error): 



Code Phase Error, d = 



D/ Areas 




2(A+A 2 ) 



[00114] 



Two situations of relevance to the region over which phase error may 



Code Phase Error, d = 



D(A 3 -A 2 -A 1 -A 4 ) - DA } 




2(A i+ A 2 ) 2(A,+A 2 ) 



(b) Delayed to Edge of Linear Region (maximum measurable phase error): 
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Code Phase Error, d = ^A 2 +A x -A>) = ~D(2A 2 + A ) = D . A 2 = D 

2{A,+A 2 ) 2{A l+ A 2 ) 2 L {A { +A 2 Y 2 V 

where a = A 2 /(A 2 + A x ) =w 2 / D and w 2 is the width of A2. 

[00115] The preceding equations for phase area can also be expressed in terms of 

the widths of each area, provided the areas take on the same height, h. Let w/, wj, W3 and 
W4 be the respective widths of Ai 5 A2, A3 and A4 as shown in FIG. 24B. Furthermore, 
each area is of identical height h. Substituting, A, = kw{, (i=l-4), D = w\ + W2 and 
substituting a ~w 2 l D ~ w 2 l(w x + w 2 ) gives the maximum and minimum phase errors 
than can be accurately measured. 

(a) Advanced to Edge of Linear Region (minimum measurable phase error): 

Code Phase Error, d - -— (1 - a) = (-w, / 2) ; 

2 

(b) Delayed to Edge of Linear Region (maximum measurable phase error): 

Code Phase Error, d = — (1 + a) = (w 2 + w, / 2) . 



[00116] In the linear region, at any given alignment of pulse signal with PN code, 

the sum of the areas is a constant. This is represented by: 



Code Phase Error, d 

2(A { +A 2 ) 2Dh 



where C is the constant, and h is the height areas Ai and A2. 
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[00117] Now assume there is a shift a small increment x to the left. An additional 

positive area of hx is "acquired" and a negative area of - hx is "lost", for a net positive 
gain of 2hx . This yields: 

D(C + 2hx) 



Code Phase Error, d - 



2Dh 



provided that the phase displacement x is in the linear region. Differentiating with 
respect to x yields the slope of the code phase error d: 

d(Code Phase Error, d)/dx = D2h 1 2Dh = 1 . 



[00118] This result shows unity slope, as is desired. Advantageously, the chosen 

phase error equation indeed measures a true code phase error, even for the more difficult 
to analyze case of asymmetric pulses. It will be appreciated that the utilization of the 
areas A] and A2 in the denominator of the above equation is of significant benefit to 
obtaining a true measure of phase error that is extremely accurate for small phase offsets. 
It may readily be appreciated that while division by other quantities, such as other areas 
and combinations of areas or signal power level determined from a prompt PN code's 
accumulation may be employed, not all such quantities provide proper scaling. Only 
areas that coincide with a level transition (or occur at the same frequency as level 
transitions, or are somehow made proportional to the number of level transitions) may be 
used. Ai and A 2 are areas that coincide with a level transition. Changes in phase error 
effect only the summations that straddle the level transitions; the summations over other 
windows (e.g., non-transition windows and windows such as A4) are relatively un- 
effected. In other words, a change in phase error will result in a change in the summation 
value in the numerator of Eq. (8) or Eq. (9) that is proportional to the number of level 
transitions as well as the size of the phase error. To achieve accurate scaling and to 
produce true phase error, the summations are divided by another summation that is also 
proportional to the number of level transitions. In addition, the summation by which you 
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divide should be designed to be independent of the exact location of the level transition 
relative to the summation window. This is why absolute values are used in the 
summation in the denominator of Eq. (8) or Eq. (9); they cause the data to take a positive 
sign on both sides of the transition so that all data is positive regardless of the transition's 
location. Finally, to complete the scaling, and to get the true phase error, a multiplication 
by a selected scaler, D/2 is employed, where D is the width of the window used in the 
summation in the denominator of Eq. (8) or Eq. (9). 

[00119] The use of Ai and A2 may have an additional advantage in that they may 

better account for the effects of signal attenuation and rounding (caused by limited 
receiver bandwidth) in the vicinity of the transition. Because the areas A\ and A 2 are 
near the transition (and, for the most part straddle the chip transition), they contain all the 
information as to the whereabouts of the transition's location (e.g., the phase error). 
Therefore, it is advantageous to divide by the absolute areas of Ai and A2 when 
normalizing the phase error to produce a true phase error (of unit slope). In other words, 
as the samples in the vicinity of the transition are summed, the areas Ai and A2 provide 
an accurate estimate of the magnitude of the samples despite the variability in the 
amplitudes and transition due to filtering effects. 

[00120] From the above analysis it is possible to plot the error curve for the 

measured phase error. The code phase error d, 322 for the asymmetric pulse signal is 
shown in FIG. 24A. Within the linear region from - (w 2 1 2) to (w 2 + w x 1 2) , the method 

accurately measures code phase error d, 322. Zero error, indicating perfect alignment of 
the pulsed (reference) signal with the received PN code, occurs when unit-sloped curve 
541 crosses the horizontal axis at 546. The pulses in FIG. 24B are shown in perfect 
alignment with the chip edge L-H transition 547. Errors to the left, denoted -x (543), 
occur when the phase of the pulsed (reference) signal leads the PN code phase. Errors to 
the right, denoted +x (544), are indicative of a lagging pulsed (reference) signal. Notice 
that when the lead is sufficiently large, the phase error curve 542 becomes zero. As will 
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be shown, this has the advantage of eliminating the effect of multipath when the pulses 
sufficiently lead the multipath signal (multipath, naturally, lags so the pulses when 
aligned with direct signal lead the multipath). It is, however, undesirable to operate the 
GPS receiver 10 such that the direct signal is in region (542) of the curve. Because the 
measured code phase error here (to the left of-W//2) is zero, there is no reliable 
indication of true phase error and loss of phase lock may occur. That is, in this region 
even if there were an error a tracking loop would not be able to correct for it. It should be 
noted that the search and acquisition logic of the GPS receiver 10 starts the initial 
tracking at a code phase error preferably close to zero, but if biased, biased to the right 
(lagging) so that proper pull-in of the track loop can occur due to an always positive 
phase error 540. 

[00121] Turning now to FIG. 25 where the code phase error curve 550 is shown 

► 

and also multipath phase error curves 551 and 552. It may readily be observed that the 
asymmetric phase error is beneficial for reducing multipath. As the figures depict, the 
multipath is always delayed with respect to the direct signal as shown by the offset 553. 
The multipath 551 in FIG. 25A is of same polarity as the direct PN signal, but it is 
equally likely that the multipath will be of opposite polarity to the direct signal as shown 
by 552 of FIG. 25B. Since multipath enters the system linearly, superposition holds. 
That is, the total phase error is the sum of code phase error 550 from the direct PN signal 
and the code phase error 551 or 552 from the multipath signal. When the multipath 
signal is delayed by more than the duration of wi, it no longer has an effect on the zero 
crossing point 554 of the direct signal's code phase curve. Since the receiver operates to 
maintain near zero code phase error, the adverse effects of the multipath are eliminated 
for multipath delayed by more than w { . 

[00122] It may further be appreciated that another anomaly arises, however, if the 

multipath is delayed by slightly more than one chip. In FIG. 24A there is shown a 
positive region 545 in the code phase error curve that is slightly more than one chip 
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advanced from the zero phase error point, 546. It may now readily be appreciated based 
upon the teachings herein, that this positive region may also be removed by adding 
additional pulses in the vicinity of the non-transition pulse A3. The goal is to make the 
composite set of pulses on a non-transition similar in shape to that of the composite set of 
pulses on a transition, but still keeping the total area of the non-transition pulses equal to 
A 3 . This may be accomplished with an additional set of positive and negative pulses, 
denoted as A5 and A 6 that have equal areas (opposite polarity) and thus cancel. Due to 
the area cancellation, all the derivations relating to FIGS. 22A-22C and FIGS. 23A-23B 
still hold. 

[00123] Turning now to Figures 26A and 26B, FIG. 26B depicts a set of pulsed 

windows for generating the asymmetric code phase error measurement in accordance 
with an exemplary embodiment. FIG. 26B shows the addition of two pulses with equal 
areas A5 and A 6 but opposite signs. As is apparent in the figure, the shape of the 
composite of pulses (Ai, A2 and A4) in the region of a transition is now similar in shape 
to the composite of pulses (A3, A5 and A6) in the region of a non- transition. FIG. 26A 
shows the resulting phase error curve. The phase error curve is virtually unchanged 
relative to FIG 24A over a region from one chip advanced 561 to one chip delayed 562. 
The left-hand portion 563 of the code phase error curve in FIG 26A now remains zero, 
even in the region 564 where the measured code phase error 545 of FIG 24 A was shown 
to be non-zero. This is a result of the additional pulses A4 and A6 that cause the non- 
transition composite pulse set to balance with the transition composite pulse set when the 
reference phase leads by a chip or more. Therefore, to reject multipath of greater than on 
chip delay and generate an asymmetric phase error measurement, a set of pulsed windows 
in accordance with an exemplary embodiment as depicted in FIG. 26B may be utilized. 

[00124] So far, the multipath mitigating phase error concept of an exemplary 

embodiment has been presented in graphical terms where pulses are depicted having 
areas Ai, A2, A3, A 4 , A 5 and A6. Although this has been beneficial for understanding the 
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concept, the actual preferred embodiment involves processing of discrete digital data. 
Recall that in an exemplary embodiment, sums over windowed sections of sampled data 
are computed rather than areas, which were employed for illustrative purposes. The 
windows are aligned in time according to the reference code phase (e.g. 324 Figure 14) 
generated by the receiver 1 0. With the disclosure and teachings presented thus far, Eq. 
(9) may now be stated with greater detail of the summation as Eq. (10) below. The code 
phase error, d, 322 is given by 



2>iu2(*)*(*)+ 2>3u5<*)«y(*)+ 2> 6 (*)*(*) 

_ r A: g 1 wtndows ^A2 MT - ntf0WK keA4^ in(jom k e A 3 ^ indow5 ^>A5 ^ n{ j OVirs windows 1 

~2 L Y\S{ki J 



Ira 4\ i i AI 

n. t /i i windows winows 



(10). 



[00125] In this instance, the notation Al v/tndows U A2 wWowj denotes the composite of 

the Ai and A2 windows and A2>^ ndov/s u A5 windows denotes the composite of the A 3 and A5 
windows. The notation k e A4 windows , for example, means that the received signal, S(k), is 
sampled at all points k such that k falls within the windows of area A4. Note that the 
Endows > A2 windows and A4 windows all occur in the vicinity of a H-L transition or L-H 

transition and the A3 windom , AS^*^ and A6 mndows occur on non-transitions as shown in 
FIG 26B. Also, as shown in FIG 26B, D is the combined width of the Ai and A2 
windows. The multipliers, a i9 (i = 1 u 2,4,3 kj 5,6) take on values of +1 or -1, and are 

defined as follows: 



a iw2 = 



+ 1 if L - H transition of reference PN code, 
- 1 if H - L transition of reference PN code; 



^4 = 



- 1 if L - H transition of reference PN code, 
+ 1 if H - L transition of reference PN code; 



-41 - 



COU-0011 




4-1 if non - transition and polarity of reference PN code is positive, 
- 1 if non - transition and polarity of reference PN code is negative; 



^6 = 



- 1 if non - transition and polarity of reference PN code is positive, 
+ 1 if non - transition and polarity of reference PN code is negative. 



[00126] 



In an exemplary embodiment, the prompt PN reference code and the 



conventional V 2 chip early l / 2 chip late reference code (E-L PN code) (See Figures 5 and 
6) are utilized as a means to compute the '+1 ' or '-1 ' sign represented by a . The E-L 
reference code is positive on a L-H transition, negative on a H-L transition, and zero on a 
non-transition of the prompt PN reference code. When E-L is zero, the prompt code 
dictates the sign of the pulses with areas A3, A5 and A6, otherwise, the E-L code dictates 
the sign of pulses with areas A t , A2 and A4. The input PN code sequence, S(k) (e.g., FIG. 
17B), is sampled in the various windows and appropriately multiplied by the correct sign 
based on the E-L and prompt reference codes. Summations are formed and the phase 
error is computed according to Eq. (10). FIG. 27 depicts a flow chart illustrating an 
exemplary process 100 for phase error computation and utilization of the phase error to 
control the reference phase and to align it with the received PN phase. 

[00127] Observation of FIG. 27 brings to bear yet another beneficial aspect of an 

exemplary embodiment of the invention. Specifically, a count 570 and 571 is maintained 
for both transitions and non- transitions. It will be appreciated from the table of FIG. 15, 
that the number of transitions (either H-L or L-H) does not necessarily balance the 
number of non-transitions for C/A code. This is phenomenon is particularly true for 
satellites with PRN IDs 7, 8, 15, 17, 21, 22, and 24. It should be further appreciated that 
the asymmetric set of pulses derived thus far for multipath mitigation presumes that the 
number of polarity transitions is equal to the number of non-transitions. However, as is 



-42- 



COU-0011 



evident from the table, they are not always equal. Moreover, it will be seen that if the 
inequality is not appropriately accounted for, a code phase error results. 

[00128] One way to overcome this limitation is disclosed in yet another exemplary 

embodiment that maintains separate accumulated sums of transitions for the areas of 
transitions and the areas of non-transitions. 



Let SumT= sum of product of PN code with a over areas Ai, A 2 and A4, and 
let SumN = sum of product of PN code with a over areas A3, A5 and A6 

[00129] At, A2, and A4 occur in the vicinity of a transition. Let Nt be the number 

of such transitions within the accumulated summation. A3, A5 and A6 occur in the 
vicinity of a non-transition. Let Nn be the number of non-transitions within the 
accumulated summation. Eq. (10) can then be modified as: 



PI 2[(Nt I Nn)SumN + SumT] 

C AU (11)- 

Sum Abs 



where 



SumT= JV^ $> 4 (*)S(*) 

keAl ^(10^ <JA2 windows k&AA windows 



SumN= 5>3w5(*)S«+ 5> 6 (*)S(*) 

ksA3 widows \JAS windows keA6 windows 



SumAbs = £|S(Jfc)| 



k€.A\ windows ^^7. winows 



-43 - 



COU-0011 



[00130] In an exemplary embodiment, the ratio Nt/Nn provides compensation to 

the summation, SumN, to account for the inequality between the transitions and non- 
transitions. Obviously, after reading this disclosure, one skilled in the art could devise 
other approaches to compensate for the lack of balance between the transitions and non- 
transitions. For example, in yet another embodiment, logic could be deployed that 
stopped summing Su m T and SumN when their respective transition, Nt, or non- transition, 
Nn 9 reached a certain value. This value could be the minimum that either Nt or Nn could 
obtain if summations were continued to the end of the accumulation interval. In yet 
another approach of an exemplary embodiment an accumulator could be employed to 
keep track of total number of sample points utilized within each summation SumT and 
SumN with the realization that the number of sample points is proportional to the number 
of transitions and non-transitions, respectively. 

[00131] Regardless, of the methodology, it will be appreciated that an exemplary 

embodiment of the invention provides compensation for unequal transition and non- 
transition densities in a pseudorandom code. Moreover, an exemplary embodiment of the 
invention removes the dependence of the code phase delay, d, on the transition type (HL 
or LH) by defining one transition window and using the S'(k) waveform to set the 
polarity of the window. Further yet, the concept may be extended to unequal HL and LH 
transitions by utilizing unique windows for the non-transition, HL transition and the LH 
transition cases and counting the corresponding occurrences with Nn 9 Nhl and Nlh, 
respectively. In one embodiment a similar shaped window (only differing in sign) for the 
HL and LH transition windows and sets Nt = Nhl + Nlh is employed. 

[00132] As a further point, although an exemplary embodiment normalizes the 

phase error measurement to achieve unit slope, it should be appreciated that this does not 
have to be the case to take advantage of the compensation approach disclosed herein. 
That is, for example, the denominator in Eq. (1 1) or even the term "Z)/2" could be 
omitted and still yield a quantity that is proportional to the code phase error, d. The 
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feedback control loop could readily be configured to use the proportional term rather than 
true error d. 

[00133] It will be appreciated that the satellite systems as discussed herein may 

include but not be limited to Wide Area Augmentation System (WAAS), Global 
Navigation Satellite System (GNSS) including GPS, GLONASS, GALILEO, and other 
satellite ranging technologies. The term WAAS here is used as a generic reference to all 
GNSS augmentation systems which, to date, include three programs: WAAS (Wide Area 
Augmentation System) in the USA, EGNOS (European Geostationary Navigation 
Overlay System) in Europe and MSAS (Multifunctional Transport Satellite Space-based 
Augmentation System) in Japan. Each of these three systems, which are all compatible, 
consists of a ground network for observing the GPS constellation, and one or more 
geostationary satellites. 

[00134] The disclosed invention can be embodied in the form of computer or 

controller implemented processes and apparatuses for practicing those processes. The 
present invention can also be embodied in the form of computer program code containing 
instructions embodied in tangible media 13, such as floppy diskettes, CD-ROMs, hard 
drives, or any other computer-readable storage medium, wherein, when the computer 
program code is loaded into and executed by a computer or controller, the computer 
becomes an apparatus for practicing the invention. The present invention may also be 
embodied in the form of computer program code as a data signal 15, for example, 
whether stored in a storage medium, loaded into and/or executed by a computer or 
controller, or transmitted over some transmission medium, such as over electrical wiring 
or cabling, through fiber optics, or via electromagnetic radiation, wherein, when the 
computer program code is loaded into and executed by a computer, the computer 
becomes an apparatus for practicing the invention. When implemented on a general- 
purpose microprocessor, the computer program code segments configure the 
microprocessor to create specific logic circuits. 
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[00135] It will be appreciated that the use of first and second or other similar 

nomenclature for denoting similar items is not intended to specify or imply any particular 
order unless otherwise stated. 

[00136] While the invention has been described with reference to an exemplary 

embodiment, it will be understood by those skilled in the art that various changes may be 
made and equivalents may be substituted for elements thereof without departing from the 
scope of the invention. In addition, many modifications may be made to adapt a 
particular situation or material to the teachings of the invention without departing from 
the essential scope thereof. Therefore, it is intended that the invention not be limited to 
the particular embodiment disclosed as the best mode contemplated for carrying out this 
invention, but that the invention will include all embodiments falling within the scope of 
the appended claims. 



) 
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